Eliashberg Analysis of Optical Spectra Reveals Strong Coupling of 
Charge Carriers to Spin Fluctuations in Superconducting Iron Pnictides 
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The temperature and frequency dependences of the optical conductivity of Co and Ni-doped 
BaFe2As2 are analyzed and the electron-boson spectral density a^F{uj) extracted using Eliashberg's 
formalism. The characteristic energy of a large peak in the spectrum around 10 meV coincide with 
the resonance peak in the spin excitation spectra, giving compelling evidence that in iron-based 
superconductors spin fluctuations strongly couple to the charge carriers and mediate superconduc- 
tivity. In addition the spectrum is found to evolve with temperature towards a less structured 
background at higher energies as in the spin susceptibility. 
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Whenever a new family of superconductors is discov- 
ered, among the first questions posed are about the mech- 
anism of superconductivity P, Q • In conventional super- 
conductors, phonons mediate the attractive interaction of 
two electrons forming Cooper pairs. For heavy fermions, 
cuprates and some organic superconductors, magnetic 
Cooper-pairing mechanisms have been proposed as a can- 
didate |3|. However, intrinsic complications have pre- 
vented the general acceptance. The situation of the novel 
class of iron based superconductors is comparably unset- 
tled Nevertheless, the proximity of superconductiv- 
ity and antiferromagnetism in the phase diagram , the 
weak electron-phonon interaction [y,01 and the resonance 
peak in the spin-excitation spectrum d, Q support the 
hypothesis of a magnetic interaction leading to supercon- 
ductivity [iMl- 



The Eliashberg spectral function a'^F{u!) quantifies the 
boson exchange effect and is a good way to discrim- 
inate between the candidate mechanisms. For half a 
century, current- voltage characteristics obtained by tun- 
neling spectroscopy arc utilized to provide detailed and 
accurate information on the phonon exchange for most 
conventional superconductors [l3 , 14 1 . In strongly corre- 
lated systems, such as high- Tc cuprates, the inversion of 
optical data is commonly used to extract the bosonic ex- 
citation spectra 15l-tl7l| . The analysis and interpretation. 



however, is not straight forward since a certain complex- 
ity in understanding the spectral signature of bosonic 
modes arises from the joint mechanisms in these partic- 
ular systems. In this regard, Carbotte and collaborators 
succeeded to modify the kernel a^F^ou) of the Eliash- 



berg theory by introducing the nearly antiferromagnetic 
Fermi-liquid model where the exchanged bosons are de- 
scribed as antiferromagntic spin fluctuations, and the re- 
sulted optical resonance tracked very well the tempera- 
ture evolution of the spin resonance seen in neutron scat- 
tering 
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In this Letter, we reported a detailed analysis of our 
optical spectra obtained on Ba(Feo.92Coo.o8)2As2 single 
crystals. The inversion of the frequency dependent op- 
tical scattering rate r~^(a;) reveals that the coupling of 
charge carriers to bosonic modes has an optimum peak 
around 10 meV, with a coupling constant A = 4.4 right 
above Tc- With increasing temperature, this peak be- 
comes broader and moves to higher energies. These 
bosonic spectral signatures closely resemble the dynami- 
cal spin susceptibility x"{^) observed by neutron scatter- 
ing, indicating a magnetic mediation mechanism in the 
novel iron-pnictide superconductors. 

We have measured the optical reflectivity of Co-doped 
BaFe2As2 single crystals over a wide frequency and tem- 
perature range as described in detail in Rcf. llSl . The sam- 
ples are well characterized and exhibit a superconducting 
transition at Tc = 25 K [l^ . Via Kramers-Kronig anal- 
ysis we calculate the complex conductivity a = ai + ia2 
which is further analyzed by the extended Drude model in 
order to obtain the frequency dependent scattering rate 
l/T{io) and mass enhancement m*(cj)/mf, = 1 + X{lu) 
compared to the bandmass mi,. The results arc plotted 
in Fig. [T] for different temperatures. 

In principle, optical data are encoded with informa- 
tion on the microscopic interaction between the charge 
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FIG. 1; (Color online) Results of the inversion calcula- 
tions of the electron-boson spectral function a^Flco) for Ba- 
(Feo.92Coo.o8)2As2 at different temperatures in the normal 
state, (a) Measured conductivity spectra (solid lines) com- 
pared with computational results (dashed lines), (b) Fre- 
quency dependent scattering rate obtained by the extended 
Drude analysis of the conductivity spectra plotted in the up- 
per panel (solid lines) compared with the calculated l/r(aj) 
according to Eq. ([T| (dashed lines) with a constant contribu- 
tion of 304 cm~^ due to impurities, (c) The corresponding 
electron boson spectral density a^F{uj) calculated by Eliash- 
berg formalism. 



carriers. For an clcctron-boson system, the Eliashberg 
equations apply and a Kubo formula can be used to get 
the infrared conductivity from the electron-boson spec- 
tral de nsity once the quasiparticle self-energy I](aj,T) is 
known [20[. The opposite direction turns out to be even 
more challenging since we have to invert an equation of 



1 1 r°° 

-— = +/ K{uj,n;T)a^F{n)dn , (1) 

T[LO) Tiinp Jo 

where l/rimp denotes a constant scattering rate due to 
impurities; the normal state kernel is given by [2l| : 



K{LU,n;T) = - 



2u! coth 



(2) 



(w - ri) coth <! ^ I + (w - il) coth 1 — — ^ 
^ ' ^2kbT/h}^ ' \2kbT/h 



Several methods have been suggested to extract the in- 
formation on the electron-boson spectral density a'^F{ui) 
from the optical scattering rate, such as singular value de- 
composition, maximum entropy method and least square 
fit; a detailed discussion of advantages and limitations of 
these numerical inversion techniques is given in Ref. . 
Here we have approached the deconvolution by an unbi- 
ased maximum entropy method, similar to jl7j |. and plot 
the outcome in Fig. [IJc) for various temperatures in the 
metallic state of Ba(Feo.92Coo.o8)2As2. 

Note, when a gap in the density of states opens be- 
low Tc, the present analysis of T^^(a;) becomes mean- 
ingless; hence we have to restrict ourselves to T > Tc- 
However, it is safe to assume that electron-boson cou- 
pling makes a strong impact on the spectra already in 
the normal state just above Tc, and it is this quantity 
which will determine Tc- Before starting the discussion, 
we would like to emphasize that the Eliashberg inversion 
applied here is based on single band systems with infinite 
band width while most of materials actually have finite 
band width. The iron pnictides, on the other hand, are 
certa inly multiband systems (see for instance Refs. (loj 
and [l^l). In the normal state the optical conductiv- 
ity of such a multi-band system is just the sum of the 
various band contributions to this conductivity, provided 
the interband optical transitions are zero or sufficiently 
small, as was suggested by van Heumen et al. for 
the Ba(Fei_a;Coa;)2As2 class of materials. As is generally 
believed, the inelastic scattering is dominated by inter- 
band transitions due to possibly spin fluctuations and, 
thus, a single form of the electron-boson spectral density 
will enter the problem except for a possible scaling fac- 
tor accounting for a different magnitude of a^F{uj) for 
transitions between different bands. This was considered 



explicitly by Benfatto et al. [2J]. Thus, the application of 
Eq. Il]) for the deconvolution of experimental data will 
yield meaningful information about the electron-boson 
interaction in such systems and will provide an average 
a^Floj) spectrum. The shape of such a spectrum will 
still provide meaningful information on the bosons re- 
sponsible for superconductivity, for example phonons or 
spin fluctuations. A knowledge of its average magnitude 
is also equally important. Finally, the use of a formula 
based on infinite band width to deconvolute optical data 
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of systems with a rather narrow band width makes it 
extremely important to check the Kramer-Kronig consis- 
tency of a'^F{uj) with optical constants to exclude any 
incorrect solutions due to finite band-width effects. In 
order to demonstrate the applicability of our analysis, 
we added the calculated cti(w) and t^^[u) to Fig. [TJa) 
and (b) as dashed lines; both cases show good agreement 
between theory and experiment. This gives confidence in 
the physical relevance of the derived spectra. 

In Fig. [TJc), one sees a clear temperature dependence 
to the recovered electron-boson spectra which is also in 
agreement with a second, biased maximum entropy inver- 
sion At temperatures just above Tc, a pronounced 
peak centered at 10 meV and a shoulder around 45 meV 
dominate the spectral weight below 80 mcV. When T in- 
creases, this peak moves to higher energies and smears 
out quickly as a shoulder. As a consequence, the mass 
renormalization factor A is reduced from 4.4 to 1.67 at 
T = 200 K. Since 10 meV seems to be a reasonable en- 
ergy for phonon excitations, at first glance, it is tempting 
to consider a phonon mechanism for superconductivity. 
However, compared to band structure calculations 0,[z[, 
the observed spectral features are quite different. More- 
over, from calculations of the electron-boson excitation 
spectra a characteristic phonon frequency oJin can be ex- 
tracted; in our case hwin = 14.2 meV. When the phonon 
mechanism is dominant in a superconducting material, 
one can estimate the coupling strength by its ratio to T,.. 
Here, we obtain kBTc/hcuin = 0.15 < 0.25, implying a 
conventional strong coupling material; it also yields A to 
be in the range 1 — 2, according to the McMillan equa- 
tion [31. Obviously, this is much too small compared 
to our experimental result (A = 4.4). Here, we would 
like to point out that the mass renormalization factor 
A is widely reported to be A=4-5 by other studies [25| . 
Thus, we expect another mechanism to play the key role 
in mediating superconductivity in these materials. 

As mentioned above, spin fluctuations seem to be the 
natural candidate for the superconducting "glue" in iron- 
based materials. Whenever a magnetic mechanism is dis- 
cussed, the main concern is whether spin- fluctuation ex- 
change provides sufficient spectral intensity to make a 
significant impact on the electronic self-energy. Very re- 
cently, Dahm et al. [i^ succeed to establish a quantita- 
tive relationship between the charge- and spin-excitation 
spectra in high-Tc cuprates, which demonstrates that 
the magnetic interaction can generate d-wave supercon- 
ducting states with transition temperatures comparable 
to the maximum Tc observed in these compounds; in 
other words, spin fluctuations do have enough strength to 
cause superconducting transitions at high-temperature. 
In Fig. [2] we display the electron-boson spectral func- 
tion a'^F{uj) derived from optical scattering rate of Ba- 
(Feo.92Coo.o8)2As2 together with the spin excitation spec- 
trum x"i^) of Ba(Feo.925Coo.o75)2As2 obtained from in- 
elastic neutron scattering experiment Both spectra 
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FIG. 2: (Color online) Comparison of the electron-boson spec- 
tral function o^F^lu) (solid lines, left axis) derived from inver- 
sion of the optical scattering rate of Ba(Feo.92Coo.o8)2As2 and 
the spin excitation spectrum x"i^) (dots, right axis) obtained 
from inelastic neutron scattering data for Q = Qafm = 
(^, |, 1) from 01 at temperatures as indicated. 



show a resonance at approximately 10 meV at low tem- 
perature, and the peak smears out when the temperature 
increases. Such an important agreement between optics 
and neutron indicates that the charge carriers in this ma- 
terial are strongly coupled to the spin excitations. 

Our results for Ba(Feo.95Nio.o5)2As2 {Tc ~ 20 K) are 
qualitatively similar, but span only a limited energy 
range and exhibit excessive noise due to the smaller crys- 
tal size [l^, [l^ . Inelastic neutron scattering indicates a 
resonance peak in the spin excitation spectrum around 
7 meV corresponding to the lower energy scale in this 
material. Yang et al. |27l | performed a similar analysis on 
K doped BaFe2As2 and found two maxima of a^F{uj) in 
the range below 30 meV. Although different in detail, the 
overall accord gives us confidence that our observations 
reveal a general behavior in this class of materials. 

In conclusion, we have analyzed the temperature and 
frequency dependences of the optical properties of doped 
BaFe2As2 via Eliashberg theory. We obtained the 
electron-boson spectral density a'^F{uj) which exhibits a 
characteristic peak around 10 mcV. This coincides with 
the resonance peak in the spin excitation spectrum and 
gives evidence that in iron-based superconductors spin 
fluctuations strongly couple to the charge carriers and 
mediate superconductivity. Also there is a strong evolu- 
tion of the spectra with temperature which agrees with 
the known spin fluctuation spectra and this would not be 
the case in a phonon mechanism for which the a'^F(uj) is 
expected to remain independent of temperature. 
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